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Summary 

An experimental and analytical study to investi- 
gate the effect of spanwise curvature on flutter was 
conducted at a Mach number of 0.7. Two series of 
rectangular planform wings of aspect ratio 1.5 and 
curvature ranging from 0 (uncurved) to 1.04 ft -1 
were flutter tested in the Langley Transonic Dynam- 
ics Tunnel. One series consisted of models with an 
NACA 65A010 airfoil section; the other, of flat-plate 
cross-section models. Flutter analyses were con- 
ducted for correlation with the experimental results 
by using structural finite element methods to perform 
vibration analyses and two aerodynamic theories to 
obtain unsteady aerodynamic load calculations. The 
experimental results showed that for one series of 
models the flutter dynamic pressure increased signif- 
icantly with curvature, whereas for the other series of 
models the flutter dynamic pressure decreased with 
curvature. The flutter analyses, which generally pre- 
dicted the experimental results, indicated that the 
difference in behavior of the two series of models 
was primarily due to differences in their structural 
properties. 

Introduction 

Wing and fin designs that are curved (out of 
plane) spanwise are being used for high-speed mis- 
siles at low altitude. These surfaces are curved span- 
wise so that in the prelaunch position they can fold 
against the cylindrical missile body (fig. 1). This 
allows the missile to fit into a launch tube. The pos- 
sibility of fin flutter is a concern because of the high 
dynamic pressures at low altitude. 

The scientific literature revealed no information 
relative to the influence of wing spanwise curva- 
ture on flutter characteristics. The objective of the 
present study was to provide analytical and exper- 
imental data on generic configurations. Two series 
of rectangular planform models of aspect ratio 1.5 
were flutter tested in the Langley Transonic Dynam- 
ics Tunnel (ref. 1). The only difference between 
models in each series was in curvature (as measured 
by the reciprocal of the radius of curvature) which 
ranged from 0 (uncurved) to 1.04 ft -1 . The two 
series of models differed in the airfoil section — the 
models in one series had NACA 65A010 airfoil sec- 
tions (fig. 2), whereas the models in the other se- 
ries had flat-plate airfoil sections with rounded edges. 
Flutter analyses were conducted for comparison with 
experimental results using both the planar lifting 
surface method termed kernel function and a non- 
planar doublet-lattice method to calculate the un- 
steady aerodynamic loads. 


Symbols 

A peak amplitude of subcritical response 
data, V 

b wing semispan, measured along the arc 

length, 2 ft 

b 0 wing half-chord, 0.666 ft 

C curvature 1/R, ft -1 

EAL engineering analysis language 

/ frequency, Hz 

GVT ground vibration tests 

g structural damping, In 

M Mach number 

M 0 model mass excluding support, slugs 

N number of cycles 

q dynamic pressure, psf 

R radius of curvature, ft 

V velocity, fps 

Vi flutter speed index, — r~~f= 

UaboVV 

x 0 reference displacement amplitude, ft 

xjy displacement amplitude after N cycles, ft 

p density, slugs/ft 3 

p mass ratio, 

nb 0 bp 

u> a first torsion mode frequency, rad/sec 

Subscripts: 

a analysis (calculated) 

m measured 

F flutter 

Apparatus 

Wind Tunnel 

The wind-tunnel tests were conducted in the 
Langley Transonic Dynamics Tunnel (TDT) (ref. 1). 
The TDT is a continuous-flow, single-return wind 
tunnel with a 16-ft square test section (with cropped 
corners) having slots in all four walls. The tunnel is 
equipped to use either air or Du Pont Freon- 12 as 
the test medium (Freon-12 was used in the present 
study) at pressures which vary from near vacuum to 
slightly below atmospheric. The flow is generated 
by a motor-driven fan. The range of Mach number is 
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from 0 to 1.2. Both the density and the Mach number 
are continuously controllable. The TDT is equipped 
with four hydraulically activated bypass valves that 
are used to rapidly reduce the Mach number and 
dynamic pressure in the test section when flutter is 
encountered. 

Models 

Geometry. Eight semispan models — four with 
NACA 65A010 airfoil sections and four with flat- 
plate cross sections — were tested. Values of curva- 
ture (the reciprocal of the radius of curvature) were 
0 (uncurved), 0.625, 0.787, and 1.04 ft -1 . The mod- 
els had a semispan or arc length of 2.0 ft and a 
chord of 1.33 ft resulting in an aspect ratio of 1.5. 
The planform geometry and the radii are shown in 
figure 3. 

Construction. The flat-plate models were con- 
structed of 0.09-in-thick 6061-T6 aluminum-alloy 
sheet. The leading and trailing edges were rounded, 
and the aluminum sheet was bent to the appropriate 
curvatures. The NACA 65A010 airfoil section models 
were constructed of 0.09-in-thick 6061-T6 aluminum- 
alloy plate to which a lightweight (1.5 lb/ft 3 ) foam 
was attached and shaped to give the desired airfoil 
shape. The foam was attached to the curved alu- 
minum in precut strips with double-back adhesive 
tape, and the surface roughness was smoothed with 
a mixture of latex paint and filler material. High- 
strength tape was added along the leading and trail- 
ing edges of the models and also on the wingtip to 
avoid foam separation due to airloads at these critical 
locations. A 0.2-in-wide strip of No. 30 Carborun- 
dum grit was added along the 5-percent chord line 
of the upper and lower wing surfaces to provide the 
desired flow transition characteristics. The grit size 
and location were based on previous experience with 
other wind-tunnel models. 

Each model was instrumented with two strain- 
gauge bridges that were attached to the aluminum 
plate near the wing root. The two bridges were 
oriented to be sensitive to bending and torsional 
strains, respectively. 

Model Mount 

A photograph showing the uncurved 65A010 
model mounted in the wind tunnel is presented in 
figure 4. The models were cantilever mounted out- 
side the tunnel-wall boundary layer on a 4-in-deep 
I-beam support fixture attached to a remotely con- 
trolled turntable. The turntable provided the capa- 
bility of changing the wing model angle of attack dur- 
ing the test (fig. 5). A 4- by 3-ft splitter plate was 


mounted to the support fixture to provide a reflec- 
tion plane. Wool tufts 3 in. long were attached to the 
splitter plate for flow visualization. During wind-on 
tests, the tuft patterns indicated there was relatively 
undisturbed flow over the splitter plate. 

Model Structural Modes 

Measured Structural Modes 

The first four natural frequencies (first bending, 
first torsion, second bending, and second torsion) 
were measured for each model. The models were ex- 
cited by using a pulsed jet air shaker, and frequencies 
were obtained by analyzing the output signals from 
the strain-gauge bridges with a frequency analyzer. 
Measured frequencies for each model are presented in 
table 1. For the measured frequencies, the frequency 
ratio (/ 2 // 1 ) decreased with curvature for each set of 
models. During the ground vibration tests (GVT), a 
change in character of the first torsion mode shape 
was noticed. Due to the inboard torsional motion 
of the wing associated with the first torsion mode, 
the outboard motion became primarily chordwise as 
model curvature increased (fig. 6). This phenomenon 
was noticed during the GVT measurements and was 
observed to be typical behavior for the curved mod- 
els. Structural damping values for the first bending 
and first torsion modes were obtained from the GVT. 
For the NACA 65A010 models these values are pre- 
sented in table 1 and were used in calculating the 
flutter results. 

Analytical Structural Modes 

A dynamic analysis using the engineering analy- 
sis language (EAL), finite element, structural analy- 
sis program (ref. 2) provided mode shapes, general- 
ized mass, and natural frequencies for the first four 
modes of vibration. The mode shapes and general- 
ized masses were used in the flutter analyses, and the 
calculated natural frequencies f a were used for cor- 
relation with the GVT-measured frequencies f m (ta- 
ble 1). The calculated node lines for the first three 
modes of the NACA 65A010 models are shown in 
figure 7 (node lines for the flat-cross-section mod- 
els were similar) and were substantiated by measure- 
ments during the GVT. Figure 8 presents the mea- 
sured and calculated node lines for NACA 65A010 
model 2 (C = 0.625 ft -1 ). In summary, the 
calculated results correlated well with the GVT 
results. The finite element models consisted of 
108 two-dimensional quadrilateral plate elements 
having both membrane and bending stiffness (fig. 9). 
The 108 element configuration consisted of 9 ele- 
ments chordwise and 12 spanwise. Only the normal 
components of modal deflections were used in the 
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flutter analysis. These were the radial components 
for the curved models. (See fig. 10.) These deflec- 
tions were also used to assess the mode shapes. 

Flutter Analysis 

Unsteady Aerodynamic Forces 

Flutter analyses were performed by using two un- 
steady aerodynamic methods. These methods con- 
sidered the wings to be thin flat plates (no airfoil 
cross section). For the present purposes one is termed 
the kernel-function method , and the other is termed 
the doublet-lattice method even though it has its own 
kernel function that relates a lifting line of doublets 
to the downwash. The kernel-function method orig- 
inated from Watkins, Woolston, and Cunningham 
(ref. 3), and its computer program is described by 
Desmarais and Bennett (ref. 4). Notably this pro- 
gram uses the planar kernel for an uncurved wing: 
thus, each span station is analyzed as if it were part 
of an uncurved wing. In contrast the doublet-lattice 
method of Giesing, Kalman, and Rodden (ref. 5) 
includes the nonplanar kernel so that the out-of- 
plane curvature of the present models can be directly 
modeled. 

Kernel- function method. The FAST com- 
puter program described in reference 4 was used 
to obtain the kernel-function unsteady aerodynamic 
forces. This program uses a surface spline (ref. 6) to 
interpolate the displacements and slopes at the down- 
wash collocation points from the calculated mode 
shapes. Figure 11 shows the 36 collocation points — 
6 points chordwise at each of 6 span stations used for 
these analyses. 

Doublet-lattice method. The doublet-lattice 
method is directly applicable to curved panels. 
From the IS AC assembly of programs (ref. 7), 
the aerostructural interface program DLIN and the 
doublet-lattice method (ref. 5) of program DLAT 
were used. Program DLIN utilizes the surface-spline 
interpolation technique (ref. 6) to obtain deflections 
and slopes on aerodynamic boxes from deflections 
at structural nodes. Program DLAT computes the 
generalized aerodynamic forces. Figure 11 shows the 
96 uniform aerodynamic boxes — 8 boxes chordwise 
at each of 12 spanwise locations. 

Flutter Computation Routines 

Program FLUTDET of the FAST routine (ref. 4) 
was used to solve the flutter eigenvalue problem. 
Unsteady aerodynamics calculated by both aerody- 
namic theories (planar kernel function and doublet 


lattice) were used by program FLUTDET (in sepa- 
rate analyses) to calculate flutter speeds at M = 0.7 
and various densities using the traditional incremen- 
tal damping approach (V-G method). From these 
calculations, matched point solutions for flutter dy- 
namic pressure and flutter speed index were ob- 
tained. (Flutter speed index, defined as U/(o; a 6 ov //7), 
is a nondimensional quantity typically used to com- 
pare parametric changes. The flutter speed index 
largely disassociates mass, stiffness, and scale effects 
from the planform or configuration being studied.) 

Wind-Tunnel Tests 

Data Analysis 

The wind-tunnel data acquisition system (ref. 8) 
was used to calculate, record, and display tunnel pa- 
rameters and model loads. The output signals from 
the model strain-gauge bridges were monitored on 
strip chart recorders. Model frequencies and peak 
response amplitudes were determined by using a fre- 
quency analyzer. The peak hold subcritical response 
method (ref. 9) was used to predict flutter onset dur- 
ing the test. This method uses data acquired by 
the frequency analyzer, which displays a continuously 
updated frequency spectrum of the response ampli- 
tude. Flutter onset is predicted by tracking the in- 
verse of the response amplitude 1/ A versus dynamic 
pressure q for a given mode. Then 1/A is extrapo- 
lated to 0 to predict the flutter q before it is reached. 

Test Procedure 

With the wind tunnel evacuated to a low den- 
sity, the tunnel velocity was slowly increased until 
a Mach number of 0.7 was reached. Subcritical re- 
sponse data were collected at this tunnel condition. 
Next, while holding the Mach number at 0.7, the 
test-section dynamic pressure was increased in incre- 
mental steps by adding Freon- 12 into the wind tun- 
nel. At selected dynamic pressures, tunnel flow con- 
ditions were held constant, and subcritical response 
data were recorded. This process was repeated un- 
til flutter was actually encountered or until sufficient 
subcritical data were obtained to predict the flut- 
ter dynamic pressure by an extrapolation technique. 
The models were tested with the angle of attack ad- 
justed to keep the lift force equal to the weight of 
the model. This corresponded to an angle-of-attack 
range of about 2° for the NACA 65A010 models and 
about 1/2° for the flat-plate models. 

Results and Discussion 

Two flutter analysis methods were used for the 
models in this investigation — a kernel-function 
method and a doublet-lattice method. These 
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methods are identified by the aerodynamic theo- 
ries used to calculate the unsteady aerodynamic 
forces. Each method used the calculated displace- 
ments perpendicular to the model in the unsteady 
aerodynamic load calculations. The analyses for the 
NACA 65A010 models utilized measured natural fre- 
quencies and structural damping from the GVT. A 
plot showing the effect of structural damping on the 
NACA 65A010 results is shown in figure 12. The 
large values of damping had a significant effect on the 
results and thus were included in the analysis. The 
flutter analyses of the flat-plate cross-section wing 
model were conducted by using measured natural fre- 
quencies from the GVT and a value of 0 for struc- 
tural damping. Damping was observed to be much 
less than that for the NACA 65A010 models and was 
not measured for these models. 

The four NACA 65A010 models and the four 
models with the flat cross section were flutter tested 
in Freon-12 at M = 0.7. 

NACA 6SA010 Models 

Analytical results. The kernel-function results 
are presented in table 2(a) and figure 13. The flutter 
speed index shows an increase as curvature increases. 
These results also show that flutter dynamic pressure 
increases with curvature. Flutter frequency ratios, 
which are also shown in the figure, indicate there 
was not a significant change in modal participation 
as curvature increases. The kernel- function method 
incorporates the mode shape and frequency effects 
on flutter due to curvature but the kernel-function 
aerodynamic theory is a planar lifting surface rou- 
tine. Therefore, the aerodynamic effects on flutter 
due to the wing curvature are not included in this 
method. The results of these analyses show the cur- 
vature effects on flutter due to the change in unsteady 
aerodynamic loads. The change in unsteady aerody- 
namic loads is associated with the changing struc- 
tural dynamic behavior (mode shapes and frequency) 
which accompany varying curvature. This change 
in structural dynamic behavior was mentioned pre- 
viously and is characterized by a change from tor- 
sional motion to in-plane or fore-and-aft motion of 
the outboard section of the wing model. The flutter 
speed index and dynamic pressure increase as this 
change occurs because less energy is obtained from 
the airstream as the torsional motion decreases in 
perpendicular amplitude. 

The doublet-lattice results are presented in ta- 
bles 2(b) and figure 13 and show an increase in flutter 
speed index and flutter dynamic pressure as curva- 
ture increases. Changes in the structural dynamic 
behavior which accompany curvature are believed to 


cause the overall increase in flutter speed. Flutter 
frequency ratios are also presented and show no sig- 
nificant change in modal participation. These anal- 
yses, like the kernel-function analyses, indicate in- 
creases in flutter speed index and flutter dynamic 
pressure, but the increase in flutter dynamic pres- 
sure is less. Because the analyses are nonplanar and 
account for the aerodynamics on the curved surface, 
this smaller increase in flutter dynamic pressure may 
be attributed to an aerodynamic effect caused by the 
model curvature. It appears that curvature causes 
the aerodynamics to decrease flutter dynamic pres- 
sure or oppose the effect of the change in character 
of the torsional mode shape. 

Experimental results. The basic experimental 
flutter results are presented in table 3. The flutter 
dynamic pressure, flutter speed index, and flutter 
frequency ratio for each NACA 65A010 model are 
presented in figure 14. The flutter frequency ratio 
increased slightly with curvature. The flutter dy- 
namic pressure and flutter speed index also increased 
as curvature increased. The change in character of 
the first torsion mode shape to in plane at the tip 
(fig. 6) is believed to cause this effect. 

Comparison of experiment and analysis. 
Basic experimental and analytical flutter results for 
the NACA 65A010 model are given in tables 2 and 
3. The NACA 65A010 wind-tunnel test results, 
the kernel-function analyses, and the doublet-lattice 
analysis all show an increase in flutter dynamic pres- 
sure as curvature increases (fig. 15). When compared 
with the NACA 65A010 experimental results, the 
kernel-function analysis flutter prediction is conser- 
vative and becomes nonconservative as the curvature 
increases. The trend to become unconservative with 
curvature is attributed to the planar nature of the 
analysis. The doublet-lattice analysis flutter predic- 
tions remain conservative throughout the range of 
investigated curvatures because the flutter dynamic 
pressure decreasing effects of curvature are included. 

Flat-Plate Cross-Section Models 

Analytical results. The kernel-function results 
are presented in table 4(a) and figure 16. These anal- 
yses predicted an increase in flutter speed index and a 
small increase in flutter dynamic pressure. Because 
the analysis is planar, the calculated increases are 
due to the changes in structural dynamic behavior 
with curvature. 

The doublet-lattice basic analytical results are 
presented in table 4(b) and figure 16. These results 
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show an increase in flutter speed index and a rela- 
tively constant flutter dynamic pressure with increas- 
ing curvature. Here the aerodynamic effect due to 
curvature appears to have offset the structural dy- 
namic effect that is evident in the kernel-function 
results. 

Experimental results . The basic experimental 
flutter results are presented in table 3(b) and fig- 
ure 17. The figure shows the variations with curva- 
ture of the flutter dynamic pressure, the flutter speed 
index, and the flutter frequency ratio. The flutter 
speed index increased with curvature as with the 
65A010 models. The flutter frequency ratio also in- 
creased slightly with curvature. The flutter dynamic 
pressure decreased with curvature between the model 
with no curvature (C = 0) and the model with the 
least amount of curvature ( C — 0.625). From this 
point on (0.625 < C < 1.04), the flutter dynamic 
pressure is nearly constant with curvature. The flat- 
plate cross-section models did exhibit the change in 
character of the first torsion mode shape to in plane 
at the tip with curvature during the GVT yet the 
flutter dynamic pressure did not increase. Analysis 
indicates that curvature causes the aerodynamics to 
decrease flutter dynamic pressure. This effect may 
have countered the effect of the change in character 
of the first torsion mode strongly enough to cause an 
overall decrease in flutter dynamic pressure. 

Comparison of experiment and analysis. 
The flat-plate cross-section experimental results 
show a decrease in flutter dynamic pressure with cur- 
vature (fig. 18) between the model with no curvature 
(C — 0) and the model with the least amount of cur- 
vature ( C — 0.625). At higher curvature, the flut- 
ter dynamic pressure remains practically constant. 
The kernel-function analysis predicts a consistent in- 
crease in flutter dynamic pressure with curvature 
that is conservative at zero curvature and becomes 
unconservative with curvature. The doublet-lattice 
analysis compares more favorably with the experi- 
ment results than the kernel-function analysis. This 
is attributed to the effects of geometry that are con- 
sidered in the analysis. The analysis shows a very 
small increase in flutter dynamic pressure with cur- 
vature between the model with no curvature (C = 0) 
and the model with the least curvature (C = 0.625). 
This trend opposes the measured trend. For higher 
curvature, the analysis predicts a slightly decreasing 
flutter dynamic pressure with curvature that com- 
pares well with the experimental results although 
somewhat unconservative. The doublet-lattice anal- 
ysis is conservative at zero curvature and becomes 
unconservative with curvature. 


Comparison of the Two Series of Models 

The experimental results for the models with 
NACA 65A010 airfoil sections showed an increase 
in flutter dynamic pressure as curvature increased. 
The models with the flat-plate airfoil sections showed 
a decrease in flutter dynamic pressure as curvature 
increased. The flutter analyses using either kernel- 
function (planar) or doublet-lattice (nonplanar) aero- 
dynamics generally predicted the flutter dynamic 
pressure for both series of models. The doublet- 
lattice results, which included the aerodynamic ef- 
fects of curvature, showed better agreement with ex- 
perimental results. Since the analyses considered the 
wings to be thin plates, they clearly indicated that 
the difference in trend for the two series of models 
was due primarily to differences in structural prop- 
erties. In other words, airfoil shape was not shown 
to be a significant factor. 

Conclusions 

An experimental and analytical study has been 
conducted in the Langley Transonic Dynamic Tunnel 
to investigate the effects of spanwise curvature on 
the flutter characteristics of a generic wing. Two 
series of rectangular planform wing models of panel 
aspect ratio 1.5 were investigated. The first series 
of wings had an NACA 65A010 airfoil section, the 
second series of wings had a flat-plate cross section. 
For each series the same four values of spanwise 
curvature were investigated. The curvature ranged 
from 0 (uncurved) to 1.04 ft -1 . The results are 
summarized as follows: 

1. The wind-tunnel experimental results for the 
models with NACA 65A010 airfoil sections showed 
an increase in flutter dynamic pressure as curvature 
increased. This was primarily due to the change in 
character of the first torsion mode. 

2. The flutter analyses for the NACA 65A010 
models using either kernel-function (planar) or 
doublet-lattice (nonplanar) aerodynamics generally 
predicted the effects of curvature. Both analyses were 
slightly conservative at zero curvature. The kernel- 
function analyses were slightly unconservative with 
curvature. The doublet-lattice analyses remained 
slightly conservative with curvature. 

3. Comparison of the kernel-function (planar) 
and doublet-lattice (nonplanar) analyses indicated 
a structural dynamic and an aerodynamic effect on 
flutter caused by curvature of the models in this 
study. The structural dynamic effect which changed 
the character of the first torsion mode shape in- 
creased flutter dynamic pressure, whereas the aero- 
dynamic effect caused by curvature of the wing ap- 
peared to reduce the flutter dynamic pressure. 
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4. The wind-tunnel experimental results for the 
models with flat-plate airfoil sections showed a de- 
crease in flutter dynamic pressure with curvature. 
These models did exhibit the structural dynamic 
change in character of the first torsion mode to 
in plane. The aerodynamics caused by curvature may 
have countered the effect of this structural dynamic 
change to cause an overall decrease in flutter dynamic 
pressure. 

5. The flutter analyses for the flat-plate mod- 
els using either kernel-function (planar) or doublet- 
lattice (nonplanar) aerodynamics predicted results 
that were slightly conservative for the uncurved 
model. Neither analysis predicted the overall de- 
crease in flutter dynamic pressure with curvature 
that was measured for the models, and both predic- 
tions were unconservative with curvature. 

6. The flutter analyses indicated that the differ- 
ence in trend for the two series of models was due 
primarily to differences in structural properties. In 
other words, airfoil shape was not shown to be a sig- 
nificant factor. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
December 20, 1988 
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Table 1. Wing Model Analytical and Measured Properties 
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Figure 2. NACA 65A010 spanwise curvature models. 
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Figure 4. Uncurved NACA 65A010 spanwise curvature model mounted in TDT. 


Turntable 



Figure 5. Sketch of model support fixture with splitter plate. 



Figure 6. Inboard torsional motion in first torsional mode (curved arrows) and associated outboard in-plane 
motion. 



(all models) 


Second bending 
(see chart below 
for location) 


Second bending node line location chart 


Mode 

Curvature, 

ft-i 

Y, ft 

1 

0.0 

1.06 

2 

0.625 

1.00 

3 

0.787 

0.97 

4 

1.04 

0.93 


Figure 7. Calculated node lines for first three modes of NACA 65A010 spanwise curvature models (Y measured 
along surface). 
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Figure 8. Experimental and calculated node lines for NACA 65A010 model. C = 0.625 ft 1 . 
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Figure 11. Aerodynamic model. 
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Figure 17. Experimental results for flat-plate cross-section models at M = 0.7. 
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